Maximum distributions of bridges of noncolliding Brownian paths 

Naoki Kobayashi,* Minami Izumi,t and Makoto Katori* 
Department of Physics, Faculty of Science and Engineering, 
Chuo University, Kasuga, Bunkyo-ku, Tokyo 112-8551, Japan 
(Dated: 27 August 2008) 

Abstract 

The one-dimensional Brownian motion starting from the origin at time t = 0, conditioned to 
return to the origin at time t = 1 and to stay positive during time interval < t < 1, is called the 
Bessel bridge with duration 1. We consider the iV-particle system of such Bessel bridges conditioned 
never to collide with each other in < t < 1, which is the continuum limit of the vicious walk model 
in watermelon configuration with a wall. Distributions of maximum-values of paths attained in the 
time interval t 6 (0,1) are studied to characterize the statistics of random patterns of the repulsive 
paths on the spatio-temporal plane. For the outermost path, the distribution function of maximum 
value is exactly determined for general N. We show that the present iV-path system of noncolliding 
Bessel bridges is realized as the positive-eigenvalue process of the 2N x 2N matrix- valued Brownian 
bridge in the symmetry class C. Using this fact computer simulations are performed and numerical 
results on the TV-dependence of the maximum-value distributions of the inner paths are reported. 
The present work demonstrates that the extreme- value problems of noncolliding paths are related 
with the random matrix theory, representation theory of symmetry, and the number theory. 
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I. INTRODUCTION 



Random walk (RW) models are important in physics, chemistry, and computer sciences. 
They can be used effectively, when we explain basic concepts of statistical physics [1], 
stochastic processes in physics and chemistry [2], and stochastic algorithms [3]. In par- 
ticular, RWs have been used to provide simple and useful models to discuss various physical 
phenomena in far-from-equilibrium, such as interface dynamics [4-6], polymer networks [7- 
9], wetting and melting transitions [10], and so on. If we take the proper spatio-temporal 
continuum limit, called the diffusion scaling limit, of the RW models, the Brownian motion 
(BM) models are obtained. By virtue of the limit procedure, the BM models are enriched 
with mathematics. The following example of a conditional BM will demonstrate this state- 
ment [11-13]. 

Consider the one-dimensional standard BM, B(t),t > 0, where (B(t)) = and 
(B(s)B(t)) = min{s,t}. The BM conditioned to stay positive B{t) > 0,0 < t < oo, is 
called the 3 -dimensional Bessel process, abbreviated as BES(3), since it is equivalent in dis- 
tribution with the radial part of the 3-dimensional BM and its transition probability density 
is given by a special case of the modified Bessel function (see Sec. II. A below and 3.3 C in 
[14], VI. 3 in [15], IV.34 in [16]). When we consider the case that it starts from the origin 
at time t = and returns to the origin at time t — 1, this conditional BES(3) is called the 
3-dimensional Bessel bridge with duration 1 starting from 0, since as illustrated by Fig.l, the 
path drawn on the (1+1) -dimensional spatio-temporal plane seems to be a bridge over the 
time axis. (Note that by the scaling property between space and time of BM, no generality is 
lost by setting the time duration be 1.) Let r(t), < t < 1, denote the 3-dimensional Bessel 
bridge. By symmetry, we can expect that the height of the bridge attains its maximum 
value with the highest probability at time t = 1/2. The probability density for r(t) G dx at 
time t — 1/2 is readily calculated as 



(see below Eq.(15)), which gives the moments ((r(l/2)) s ) = 2 s (s + 1)!! if s is even and 
((r(l/2)) s ) = {2/v/7r}2- s / 2 ((s + l)/2)! if s is odd. The shape of the present bridge is, 
however, randomly deformed, and as we can see in Fig.l the time < r < 1, at which r(t) 
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FIG. 1: A sample of path of the 3-dimensional Bessel bridge with duration 1. In this example the 
time r, at which the height of the bridge attains its maximum, is less than 1/2. 

attains its maximum, will fluctuate around the mean time (r) = 1/2. We define 



(i) 



max r(t) = r(r) 

0<t<l 



We can show that the probability density for G dh is given as 

oo 

p H(1) (h)=8j2 e" 2n2/l2 (4n 4 /i 3 - 3n 2 h), < h< oo, 

71=1 

which gives the moments 

((#i (1) ) s > = ^=^r(*/2)C(a), 8 = 0,1,2,... 
with the gamma function T(z) = J °° e~ u u z ~ l du. Here £(s) is Riemann's zeta function 



(2) 



OO 



n=l 



which is an important special function in number theory [12, 17]. In the present paper, 
as multivariate generalization of the Bessel bridge [18], we study the iV-path systems of 
the 3-dimensional Bessel bridges with duration 1, conditioned never to collide with each 
other in < t < 1; r^ N \t) = (r[ N \t), r { 2 N \t), . . . , rjf \t)), < t < 1, with the conditions 
r (A0( ) = r W(i) = and < r?°(t) < r^(*) < • • • < r^\t),0 < t < I. 

The systems of RWs with nonintersecting condition were introduced by M. E. Fisher as 
mathematical models to describe the wetting and the melting transitions and named vicious 
walk models [10]. They have been used not only to discuss dynamics of domain walls and 



melting of commensurate surfaces [19], but also to study polymer networks [8, 9], the related 
enumerative combinatorial problems [20-22], and nonequilibrium critical phenomena [23]. 
In general the noncolliding diffusion particle systems are obtained as the diffusion scaling 
limits of the vicious RW models [24-26]. In particular, the version of the vicious RW model, 
whose continuum limit is the noncolliding Bessel bridges, r^(t),0 < t < 1, is called the 
N -watermelons with a wall [9, 21, 22, 26, 27]. 

In the discrete mathematics the one- dimensional RWs conditioned to visit only nonneg- 
ative sites N = {0, 1, 2, ... } is called the Dyck paths and the asymptotics of the average 
height of the Dyck paths in the long-step limit was studied by de Bruijn, Knuth and Rice 
[28]. Recently Fulmek generalized this classical result by evaluating the asymptotics of the 
average height of the 2- watermelons with a wall [29]. In this calculation, he showed the 
fact that the number-theoretical special functions, such as Jacobi's theta function and the 
double Dirichlet series are useful to describe the asymptotics. Motivated by this important 
observation, the present authors [13] studied the N = 2 case of the noncolliding Bessel 
bridges, r^ 2 \t) = (r^\t) , r^\t)) , which is the continuum limit of the 2-watermelons with 
a wall, and clarified that this phenomenon is indeed the generalization of the relationship 
between the maximum- value distribution of the 3-dimensional Bessel bridge and Riemann's 
zeta function mentioned above. 

We will report in this paper both of the exact results and the numerical results on the 
maximum-value distributions of iV paths in the noncolliding Bessel bridges. Main exact 
result is the following determinantal expression for the maximum-value distribution of the 
outermost path (i.e. the height of the continuum limit of watermelons), 

H { N N) = max rf\t) (3) 



for the iV-path system r^'fi) = (r[ N) (t) , r { 2 N) (t) , . . . ,rf\t)); 



P(H ( ^ ) <h)= - \ } det 

2 II =i(2j - 1)! i<i.*<* 



^2 H 2 ( j+k -i)(V2nh)e~ 



2n 2 h 2 



n=—oc 



(4) 



where H n (x) denotes the n-th Hermite polynomial defined by 

I n / 2 1 /_i \k(n \n-2k 

H ( x) = n \ V { } { } 



(5) 



with [a] = the largest integer that is not greater than a. (We can see that, since H 2 (x) = 
Ax 2 -2, Eq.(4) with N = 1 and the relation p H w (h) = dP(H[ 1] < h)/dh give the result (2).) 
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FIG. 2: A sample of paths of N = 10 noncolliding Bessel bridges with duration 1, all starting from 
and returning to 0, realized by the eigenvalue process. 

The long-step asymptotics of all moments of the height of the N- watermelons with a wall 
have been fully studied for arbitrary iV > 1 by Feierl [30, 31]. We will show that our result (4) 
for the distribution functions of the continuous model is consistent with the results by Feierl 
for the moments of his discrete model. Quite recently Schehr et al. [32] showed their study 
on the same problem and others by the path-integral technique, a different method both 
from ours and Feierl's. They also reported an expression for the maximum-value distribution 
(Eq.(5) in their paper [32]), which is different from (4). We will show that the equivalence 
of these two expressions is guaranteed by the functional equation satisfied by Jacobi's theta 
function #3(2;, y). 

Dyson introduced a matrix-valued BM, M(t) = (Mj k (t)),t > 0, in the space of Hermitian 
matrices. If the matrix size is N, the N diagonal elements, Mjj(t), 1 < j ' < N, the N(N — 
l)/2 real parts of the upper-triangle elements, Ke(Mjk(t)), 1 < j < k < N, and the N(N — 
l)/2 imaginary parts of them, lm(Mjk(t)), 1 < j < k < N, are given by independent one- 
dimensional standard BMs, the total number of which is iV 2 . Dyson showed that the N 
eigenvalues of M(t) behaves as an interacting diffusion particle system on the real axis R, in 
which the long-ranged repulsive forces work between any pair of particles with the strength 
proportional to the inverse of the distance of two particles [33]. This eigenvalue process is 
called Dyson's BM model and it has been proved to be equivalent in distribution with the 
system of iV-BMs conditioned never to collide with each other [34, 35]. The correspondence 
between eigenvalue processes of matrix-valued diffusion processes and noncolliding particle 
systems has been studied [35-39]. In the present paper we will use the fact that the iV- 
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noncolliding Bessel bridges can be realized as the positive-eigenvalue process of the 2N x 2N 
matrix-valued Brownian bridge, whose distribution at each time < t < 1 is related with 
the random matrix theory [40] with a special symmetry called the class C in [41, 42]. Figure 
2 shows a sample of paths of the N = 10 noncolliding Bessel bridges realized by this 
eigenvalue process. One can imagine that it is very hard to simulate such paths all starting 
from the origin and returning to the origin with noncolliding condition by direct computer 
simulation. The present paper will demonstrate that the relationship between the random 
matrix ensembles and the noncolliding particle systems [39] provides a practical method to 
study such conditional processes effectively by computer simulations. We will report the 
numerical evaluation of the maximum-value distributions of not only the outermost path 
rpp(t), but also inner paths r^\t), k = 1, 2, . . . , N — 1. 

Both from the view-points of statistical physics and of random matrix theory [18, 40], 
the study of iV — > oo limit is interesting and important for noncolliding paths [43-45]. For 
the average value of the maximum values of the outermost path, we can read the following 
behavior from the numerical work by Bonichon and Mosbah for the watermelons with a wall 
[46], 

(H { £ } ) ~ Vl.67^-0.06, N -> oo. (6) 

Recently, Schehr et al. [32] gave an argument that the numerical data of Bonichon and 
Mosbah (6) shows only a pre-asymptotic behavior in large N and the true asymptotics 
should be 

(H ( N N) )~V2N, N^oo. (7) 

In the present paper, we report the numerical study of the iV-dependence of the maximum- 
value distributions systematically not only for the outermost path but also for all inner paths 
and discuss the N — > oo asymptotics based on our numerical data. 

The paper is organized as follows. In Sec. II. A, after giving brief explanations of the 
3-dimensional Bessel bridge and the Karlin-McGregor formula [47], we define the N- 
noncolliding Bessel bridges by giving the transition probability densities. A matrix rep- 
resentation of the symmetry class C is shown in Sec. II. B and the matrix- valued Brownian 
bridge in the symmetry class C is introduced. The equivalence in distribution between its 
positive-eigenvalue process and the noncolliding Bessel bridges is then stated. The problems 
studied in this paper is announced in Sec. II. C. In Sec. III. A the exact expression of the distri- 
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bution function of the maximum value for the outermost path, which is the same as Eq.(5) 
in [32], is derived by our method (Lemma 1). This exact expression is then transformed 
into two kinds of determinantal expressions (Proposition 2 and Theorem 4) in Sec. III. B, one 
of which is Eq.(7) given above. The key lemma in the transformation is a set of equali- 
ties between infinite series involving the Hermite polynomials (Lemma 3) derived from the 
functional equation of Jacobi's theta function $3(2;, y). The numerical study is reported in 
Sec. IV. Concluding remarks are given in Sec.V. Appendices are used to give proofs of some 
formulas. 



II. MODELS AND PROBLEMS 



A. Transition probability density of noncolliding Bessel bridges 

Let B(t),t > 0, be the one-dimensional standard BM starting from the origin; B(0) = 0. 
For any series of times t = < t\ < t 2 < ■ • ■ < tM, M — 1, 2, . . . , the probability that the 
BM stays in interval [a m , b m ] at each time t m , m = 1, 2, . . . , M, is given by 



p(B(t m ) G [a m ,b m ],m = 1,2,...,m) 

pbi rb-2 i>b M ^-t- 1 

/ dxi / dx 2 --- dx M Y[p(t m+1 -t m ,x m+1 \x m ), 

J ai Jan J a m ™_n 



'01 Ja 2 Ja M m=0 

where the transition probability density p(t, y\x) from the position x to y during time period 
t is given by the probability density of the normal distribution with mean and variance t, 

1 ( (y- X )2) 

p{t,y\x) = -^=exp|- v * - j, t>0,x,yeR. 

We consider the situation that an absorbing wall is set at the origin and BM is absorbed 
if it arrives at the wall. By the reflection principle of BM [14], the transition probability 
density of such an absorbing BM is given by 

Pebe(t,y\x) = p(t,y\x) -p(t,-y\x) 

= _J_/ e -(y-*) 2 /(2t) _ e -(»+x) 2 /(2t) I /gN 

\/27rU J 
for x, y > 0,t > 0. The survival probability of the absorbing BM starting from x > for 
the time period T is then given by 



POO 

Af(T,x)= p ahs (T,y\x)dy, 
Jo 



whose asymptotics in x/VT — > is easily evaluated as 

The transition probability density of the BM under the condition that it stays forever in the 
positive region R + = {x e R : x > 0} is then given by 

PBES(3)(*,y|x) = hm 

= -Pabs(^yk) (9) 

for rr > 0, y > 0, < t < oo. The diffusion process whose transition probability den- 
sity is given by (9) is called the 3-dimensional Bessel process, abbreviated as BES(3), by 
the following reasons. Consider the ci-dimensional BM, B^(t) = (B^t), B 2 (t), . . . , B d (t)), 
whose coordinates are given by independent one-dimensional standard BMs {Bj{t)}j =l . The 
distance from the origin of the ci-dimensional BM 

R( d )(t) = \B^(t)\ 



= y/(B 1 (t))* + (B 2 (t))* + -- + (B d (t)y 

can be regarded as a diffusion process in R + U {0} and its transition probability density is 
given by 

PbesmC*, Vis) = ?^\e-Wy™I v (f ) (10) 
for x > 0, y > 0, t > with 

d-2 

2 ' 

where I u (z) = ^^ (z/2) 2n+; y{r(n-|-l)r(z/-|-n-|-l)} is the modified Bessel function with the 
gamma function T(z) = J °° e~ u u z ~ 1 du. The process (t) is called the ci-dimensional Bessel 
process, BES(rf), [14-16]. The transition probability density (9) of the BM conditioned to 
stay positive is equal to (10) with d — 3, i.e., v = 1/2, since h/2(z) = (e z — e~ z ) / \f2nz. 
Consider the space of all configurations of iV particles in R + with a fixed order of positions, 

= {x = (xi, x 2 , ■ ■ ■ , x N ) e R+ : Xi < x 2 < ■ ■ ■ < x N }, 

which is called the Weyl chamber of type Cjv in the representation theory [48]. For x = 
(xi,X2, • • • , xn), y = (yi, Vii ■ ■ ■ , J/jv) £ W^-, t > 0, consider the following determinant 



det 

l<j,k<N 



N 

PBES(3)(t,yj\xk)] =T[ — ,A?l„ Pabs(t,yj\x k ) 
j=l J 



where the equality is given by the relation (9) and multilinearity of determinant. By the 
theory of Karlin and McGregor [47] (see also [49] with [50] and [51]), the probability that N 
particle system of BES(3)'s starting from the configuration x G can keep the order of 
N particle positions by avoiding any collision of particles for time period T > is given by 



N 



pco p 

(T,x)= / d yi ... / 
Jo Jo 



oo poo 

dyi--- / dy N det 

l<j,k<N 



Pbes(3) (T,yj\x k ) 



By the Markov property of diffusion processes, if we assume that the configuration at time 
t — 1 is fixed to be y G W^, for < t± < t 2 < 1, the transition probability density from 
x ( 1 ) g at time t\ to x^ G W% at time t 2 is given as 



pW (tl>x (i) ;t2>x (2); 



det 

l<j,k<N 



pbes(3)(i -t2,y j \x { k ) ) 



det 

l<j,k<N 



PBES^-^fl^) 



det 

l<j,k<N 



Pbes(3)(1 -ti,yj\xf ] ) 



det 

l<j,k<N 



Pabs(l - t 2 ,yj\xf ] ) 



det 

J l<j,k<N 



P^-^xflx^) 



det 

l<j,k<N 



Pabs(l -ti,yj\x^) 



Let |x| 2 = J2f=i x ] an( i define 



^(x)= n < 

l<j<k<N 



N 

x 2 j) Yl %e 



e=i 



for x = (xi, X2, ■ ■ ■ , xn) G W^t. Since we have known the asymptotics 



det 

l<j,k<N 



Pzbs{t,yj\x k ) 



r V(2iV+l)/2 e -|x| 2 /(2t) 



2^- 1 )/ 2 nf=i{r(j)r(j + i/2)} 



^(x)^(y) 



in |y| — > (the v = 1/2 case of Eq.(33) in [39]), (11) gives the following, 

pWfa.xW^x' 2 ))^ lim jf^.xWjfe.xP)) 

|y|— o y 



x det 

l<j,k<N 



exp 



,(2) 12 



+ 



X 



(1)]2 



2(1 -t 2 ) 2(1 -t x ) 



for x«, x( 2 ) G and < h < t 2 < 1, and 



p (A °(0,0;t,x) = lim p (7V) (0,x (1) ;t,x) 

|x«|-»0 

"(7r/2)^nL(2j-l)!^ MJ J P l 2t(l-t) 



(11) 



(12) 



(13) 



(14) 



for x G and < t < 1. The iV-particle system of noncolliding 3- dimensional Bessel 
bridges with duration 1 all starting from the origin is defined as the diffusion process such 
that its transition probability density is given by (13) and (14), and denoted by r^ N \t) = 
(r[ N \t) , r2 N \t) , . . . , r^\t)), < t < 1. That is, for any sequence of times t = < ti < t 2 < 
■ ■ ■ < t M < 1, M = 1, 2, . . . , and for any sequence of regions A m G W^-, m — 1,2, ... , M, 

P(r (A °(U G A m ,m = l,2,...,M) 

„ „ M-l 

= / dxW.../ rfxWnP (JV) (^,x( m );Wi,x( m+1 )). (15) 

Note that no generality is lost by setting the time duration be 1, by the scaling property 
between space and time of the present iV-particle system inherited from BM via BES(3). 
From now on we call r^ N \t),0 < t < 1, simply the N -noncolliding Bessel bridges for short. 
Remark that, if we set N = 1 and t = 1/2 in (14), hf(x) = x and Eq.(l) is obtained. 



B. Matrix-valued Brownian bridge in symmetry class C 

For N > 1, let I N be the N x N unit matrix and define the 2N x 2N matrix 

J = 




where 0^ denotes the N x N zero- matrix. Let 7i{N) and S(N; C) be the collections of all 
N x N Hermitian matrices and all N x N complex symmetric matrices, respectively. Then 
consider the following set of 2N x 2N Hermitian matrices, 

H C (2N)= \ c=\ H S ) ;H e H(N), S G S(N; C) > , (16) 

where l H denotes the transpose of H and = l S denotes the Hermitian conjugate of S. 
We can see that any element C G H C (2N) satisfies the relation 

t CJ+JC = 0, (17) 

which means that C G H C (2N) C H(2N) satisfies the symplectic Lie algebra; sym- 
bolically written as H C (2N) = sp(2N,C) n H(2N) [48]. Due to the additional sym- 
metry (17), the 2N real eigenvalues of C G H C (2N) are given in the form A = 

10 



(Ai, A2, . . . , Xn, — Ai, — A2, . . . , — Ajv), where Aj > 0, 1 < j < N. Altland and Zirnbauer 
studied H C (2N) as the set of the Hamiltonians in the Bogoliubov-de Gennes formalism for 
the BCS mean-field theory of superconductivity, with regarding the pairing of positive and 
negative eigenvalues (Xj, 1 < j < N, as the particle-hole symmetry in the Bogoliubov- 
de Gennes theory [41, 42]. They called 7i c (2N) (a representation of) the symmetry class 
C, since sp(2iV, C) is denoted by Cat in Cartan's notation [52]. 

The Brownian bridge with duration 1 starting from the origin is defined as the one- 
dimensional standard BM starting from conditioned to return to at time t — 1 and 
denoted by b(t), < t < 1. The transition probability density of b(t) is given by 

P(l - t,0\y)p(t - s,y\x) 



p B b(s,x; t,y) = 



(18) 



p(l — s, 0\x) 

for < s < t < 1, x, y e R. Let b p jk (t), < p < 2, 1 < j < k < N, and b%(t), 1 < j < k < N, 
be independent Brownian bridges with duration 1 starting from the origin. Put 



for < p < 2, and 



b%(t)/V2 i£j<k 
b%(t) if./ k 
^(t)/v/2 if j>* 

b%(t)/V2 ifj<k 
if j = k 
-b%(t)/V2 i£j>k, 



and consider the N x N matrices S p (t) = [S^ k (t)j ^ ^,0 < p < 2, and 



(«?*(*)) 



l<j,k<N 



. Then the 27V" x 2iV matrix-valued BM is defined by 



cW(*) 



< t < 1. 



(19) 



S°(t)+iA°(t) SW+iS 2 ® 
S^-iS 2 ^) -(S°(t)-iA°(t)) 

In order to define CW(t), we have used N(N + l)/2 x 3 + N(N - l)/2 = iV(2N + 1) 
independent Brownian bridges. By definition (19), C^ N \t) E 7i c (2iV),0 < t < 1, and 
C^(0) = C^(l) = 02at. That is, C^ N \t) can be regarded as a Brownian bridge in the 
N(2N + l)-dimensional space H C (2N). 

At each time < t < 1, C^fi) can be diagonalized by a unitary-symplectic matrix and 
we can obtain the eigenvalue process \^ N \t) = (X[ N \t), . . . , X^\t), —X[ N \t), . . . , —X^\t)) 
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FIG. 3: A sample of paths of the eigenvalue process '(t),0 < t < 1, of the matrix-valued 
Brownian bridge C^ N \t),0 < i, < 1, in the symmetry class C for iV = 10. 

with < \[ N \t) < ■ ■ ■ < \pp(t). Using the generalized Bru's theorem given in [39], we 
can determine the transition probability density for the positive part of eigenvalue process 
A?°(f) = (A^f), \ { 2 N \t), A#°(t)). The result is exactly the same as Eqs.(13) and (14). 
It implies that \+^(t) G W^,0 < t < 1, with probability one, and the present positive- 
eigenvalue process X^\t) is equivalent in distribution with the noncolliding Bessel bridges 
r^(t). Figure 3 shows a sample of paths of the eigenvalue process A (A °(t),0 < t < 1, for 
N — 10 generated by computer (see Sec. IV for detail). There we can see ten positive paths 
A+^(t) and their mirror images with respect to the line A = 0. The sample of paths of the 
noncolliding Bessel bridges, r^(t),0 < t < 1, shown in Fig. 2 for iV = 10 is just obtained 
by the upper half of this figure. 



C. Problems 

For the iV-path system of noncolliding Bessel bridges, r(t),0 < t < 1, we study the 
maximum values for each path attained in the time interval t e (0, 1), 

Hi N) = maxr ( k N \t), l<k<N. (20) 

The problem considered in the present paper is to clarify the statistical property of the 
random variables H^, 1 < k < N. We will report the exact expressions of the distribution 
function for the outermost path for general N > 1 and the numerical results for inner 



12 



paths. The dependence of N is studied and the asymptotics in N — > oo will be discussed. 



III. EXACT RESULTS FOR THE OUTERMOST PATHS 
A. Distribution function of 

In this section we derive an exact expression for the distribution function of the maximum 
value of the outermost path, ~P(H^ < h). In order to that first we consider the absorbing 
BM in an interval (0, h) for h > 0, in which two absorbing walls are put at the origin and 
at the position x = h. The transition probability density p^ hs (t,y\x) for t > 0, < x, y < h 
is the solution of the diffusion equation du(t,y)/dt = (l/2)d 2 u(t 1 y) / dy 2 with the initial 
condition lim t ^, u(t, y) = 5(y — x) and with the Dirichlet boundary conditions u(t, 0) = 
u(t,h) = 0,t > 0. By the method of separation of variables and the Fourier analysis, the 
unique solution is determined as (see, for example, [53]), 

h r ix 2 >^ / n 2 7T 2 \ . /nix \ . (mx \ 
Pabs(^yk) = ~ h ^exp l-—t sm (j-y) sin (j-x) (21) 

n=l v 7 

for t > 0, < x, y < h. Note that it is the different expression of the function 

oo 

y\ x ) = Yl {p(t,y\x + 2hn) -p(t,y\ -x + 2/m)} 



n=— oo 



1 oo 
n=— oo 



exp s| -^(Z/ - (x + 2hn)) 2 I - exp i ~^(y ~ (~x + 2hn)f 



which was used in our previous paper [13]. 

Consider the following two determinantal functions, 



q V0, 



) (t,y|x)= det \p ahs (t,y,\x k )}, t>0,x,yGW^ (22) 
and 

^(^ylx) = i det \p h ahs (t, yj \x k )\ : t > 0, x, y G W%, (23) 

l<3,k<N L 

where Wj^ = {x G (0, h) N : X\ < x 2 < ■ ■ ■ < x^}. By the theory of Karlin and McGregor 
[47], q( N )(t,y\x) (resp. q^\t, y|x)) is the probability for the iV-dimensional absorbing BM 
starting from x G (resp. x G W^) at time t = to survive during time period t 

by avoiding any hitting with the absorbing boundaries of (resp. W^), and to arrive 
at y G (resp. y G W^) at time t. Note that if we think that the iV-dimensional 
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vector x = (xi,X2, ■ ■ ■ ,xn) represents the positions of N particles on R + , a hitting with 
the boundary of ~W% means a hitting of the innermost particle X\ with the origin or any 
collision between neighboring particles Xj = Xj + ±, 1 < j < N — 1. Similarly a hitting with 
the boundary of means X\ = 0, or any collision of particles, or a hitting of the outer 
most particle xn with the wall at x = h. 

Consider the iV-particle system of noncolliding BES(3), R (N) (t) = (R,[ N \t), ■■■ , B$\t)) 
starting from the configuration x G at time t — 0; R^(0) = x, and arriving at the 
configuration y G at time t — 1; R^'(l) = y. Let = max 0<t<1 R^\t). Then we 

can say that P(#j£° <h) = gJ V) (l,y|x)/g( JV )(l,y|x). By definition of the iV-noncolliding 
Bessel bridges, rW(f),0 <t<l, given in Sec. II. A, we can conclude that 

P(H™<h) = hm Ct 1 ^ - (24) 
1 N ' |xHo,|yHogW(l,y|x) 1 ; 

As shown in Appendix A, the method of the Schur function expansion [39] gives the following 
asymptotics for q^ N \l,y\x) and gj^(l,y|x) in |x| — > 0, |y| — > 0, 

/9 n N/2 N 

? {A °(i,y|x) = ^-J T)!^ (x) ^ (y)x { 1 + 0(|x| ' |y|) }' (25) 



2\ JV /7T\2JV 2 " 



x £ ^p{-^|n| 2 }{^w} 2 x{i + o(M.M)}- (26) 

l<rii<n2<---<njv 

where /i^ was defined by (12). Then Eq.(24) gives the following result. 
Lemma 1. For N > 1, h > 0, 

P(lfW < fc) = c N h~ N ^ {^(n)} 2 exp[-|^|n| 2 }, (27) 

l<ni<n2<---<njv 

where C7V = 2^ 2 vr 7V ( 4W + 1 )/ 2 / nL(2j - !)!■ 



Remark 1. This expression is exactly the same as Eq.(5) in [32], which was derived by 
the path-integral method using a Selberg's integral. Here we would like to put emphasis 
on the resemblance between the summand of (27) and (14). As mentioned in Sec. II. B, 
Eq.(14) is the same as the probability density of the eigenvalue-distribution of random 
matrices in the class C (with variance t(l —t)). The exponent of the factor h- N (' 2N + 1 ) i n 
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(27) is the dimension of the space H C (2N). Another evidence to show the hidden symmetry 
of the present maximum-value problem is the following. The character of the irreducible 
representation corresponding to a partition \i of the symplectic Lie algebra is given by [48] 



det 

/ \ l<j,k<N 



' u k +N-k+l _ -(nk+N-k+l)] 



\ 



det 


~ x N-k+l _ x -(N-k+l)^ 


l<j,k<N 





(28) 



If we set rij = /ijv-j+i + j, 1 < j < N, and xj = 1, 1 < j < N, (28) gives (see, for example, 
Eq. (3.33) in [25] and Eq. (3.10) in [53]), 

sp (1 1 ••• 1) - ^ (n) 

The above observations imply that the maximum-value problems of noncolliding diffusion 
problems will be related with some enumerative problems of combinatorics in the ensembles 
of irreducible representations of symmetry, which can be regarded as discrete version of 
random matrix ensembles. The wall restriction for paths (i.e., the stay-positive condition 
for particles) is mapped into the symplectic (the class C) structure in the present case. As 
pointed out in [32], the trivial fact lim^oo P(H^ < h) = 1 for (27) gives a version of 
Selberg integral [54, 55], 

oo oo 2 N 

s e ••• e n {oKf-oM 2 } n{(M 2 e- w/2 <5} 

n\=— oo njv=— oo l<j<k<N £=1 

N 



/oo rco ly 

•••/ II {^-xirHlxje-^dxe} 
oo J-co x^akat e=1 

N 

= (2n) N / 2 N\H(2j-l)\, (29) 
which is the special case with 7 = 1 and a = 3/2 of Eq. (17.6.6) in [40]. 

B. Determinantal expressions and Jacobi's theta function 

From (27), we can obtain the following determinantal expression for the distribution 
function of 
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Proposition 2. For N >l,h> 0, 

E 



P(H { N N> < ft) = c N h- N t 2N * 1> det 

l<j,k<N 



oo 

n 2(j+fe-l) e -7r 2 n 2 /(2? l 2 ) 



n=l 



-^4 h~ N ^r N ( £-) , (30) 



where 



r N (u) = det 

l<j,k<N 



Qj+k-l 



with 0(u) = e ~™ 2u - 



The proof is given in Appendix B. 

Now we consider a version of Jacobi's theta function 



oo 

,n 2 Jin 

Z 



n=— oo 
oo 

n=— oo 

where we have set z = e" ! and q = e ym . The following functional equation is satisfied (see 
Sec.10.12 in [56], Sec.A.3.1 in [29]), 

Mw) = Mx/v, -i/y)e-™ 2/v ^-. (32) 

From this equation, we will obtain the following equalities. 
Lemma 3. For a = 0, 1, 2, • • • , 

oo , a a _j_i oo 

n « e -W+^ = £ ff a (V5Fn^ + v^C/^e-^^^ 2 , (33) 

n=— oo n=— oo 

where H a (x) is the a-th Hermite polynomial defined by (5). 
Proof. Setting x = i£ + rfu/ (2n),y = irf, £, 77 G R, in (32) gives 

By definition (31), this implies 

00 00 

^ e -7rn 2 /7? 2 +(27r^/?7 2 H-M)n _ ^ ^ e -(^r ] n+^^/ri-ii 1 u/(2 v ^)) 2 (34) 
?i=— 00 n=— 00 

We note that the generating function of Hermite polynomial is given as 



J2ux—u 

= ; n „[.c) 

nl 



n=0 
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Then (34) becomes 



oo oo 

EE 

a=0 n=—oo 



[nu) c 



-txn 2 / if +2iri£n / rf 



oo oo / . \ a 

=»e e y^- + ^ M '- 

a=0 n=—oo \ V / 

Since this equality holds for any value of u, Lemma 3 is proved. 

By setting £ = 0, 77 = h^/ljn and a = 2£ in Eq.(33) of Lemma 3, we have the equalities 



ra=l 



(35) 



for £ = 0, 1,2, Then Eq.(30) is transformed into the following other determinantal 
expression, which was announced in Introduction. 
Theorem 4. For N > 1, h > 0, 

(-1) N 



P(HjP <h) = 



det 



2 N2 1] =i(2j - 1)! i<i.*<* 



^2 H 2 (j+k-i)(V2nh)e 



2n 2 h 2 



(36) 



As special cases of (36) with N = 1 and N = 2, we have 

oo 

P(tf{ 1} < /i) = ~ £ iJ 2 (v / 2n/ i ) e - 2/l2 ™ 2 



and 



oo oo 



2 4 x 3! 



H 4 {V2n 2 h) H 6 (V2n 2 h) 



ni=— oo n,2=— oo 

Since # 2 (x) = 4x 2 - 2, # 4 (z) = 16x 4 - 48x 2 + 12 and H 6 (x) = 64a; 6 - 480x 4 + 720x 2 - 120, 
they give 

oo 

P(^ 1} < h) = (1 - 4h 2 n 2 )e- 2h2n2 (37) 

n=— oo 

and 



oo oo , 

P(^ 2) < ft) = Yl E e- 2A2 (" 2+ "2) J 1 _ l 6 /,2 n 2 + 24 ^4 n 4 + 24/l 4 n 2 n 2 _ ^6^6 

n i — — rvi n o — — rvi ^ 



ni = — OO Tl2 = — OO 



128. 



(38) 



17 



From (37) we will obtain (2). Eq.(38) is exactly the same as the result reported as Lemma 
3.1 in our previous paper [13]. 

Remark 2. Since the derivative of the distribution function P(H^ < h) with respect 
to h gives the probability density for e dh, the s-th moment of H^p, s = 1,2, ... , is 
calculated as 



where the integral by part was done. If we insert the expression (30) into (39), we have a 
determinantal expression for the s-th moment. It is essentially the same as the expression 
of Feierl given to the dominant term in the long-step asymptotics of the moment of the 
height distribution of watermelons with a wall. See the function given in Theorem 1 
of [31]. Moreover, the determinantal expression for the central limit theorem of the height 
distribution of the watermelons with a wall given by Feierl (Eq.(28) in Theorem 2) [31] can 
be identified with our second determinantal expression (36), since the functions (f>k(z) used 
there are nothing but the Hermite polynomials. We have enjoyed the perfect coincidence of 
the results obtained by the two different routes to the problem. See [13] for more detailed 
discussion on the relationship between our treatment and that by Fulmek [29] and Feierl 
[30, 31, 57]. 

IV. NUMERICAL STUDY 

A. Bessel bridge realized by eigenvalue process 

We have prepared a computer program to generate samples of paths of b(t),0 < t < 1; 
the Brownian bridge with duration 1 starting from the origin, in which each sample path is 
approximated by random walk with 10000 steps. 

Assume that we have generated three independent Brownian bridges, bj(t),0 < t < 
1, j = 1, 2, 3, by this computer program. Then we consider the following 2x2 matrix- valued 
process, 




(39) 




hit) 
h ~ ih(t) 




(40) 
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-0.5 0.5 1 1.5 2 2.5 3 

FIG. 4: Comparison of the distribution of the maximum values of the 3-dimensional Bessel 
bridge between the numerical simulation (plotted by circles) and the theoretical values (given by 
curve) . 

It is easy to see that the eigenvalue process of C^(t) is given by the positive and negative 
pair of eigenvalues A (2) (t) = (X^\t), -A? \t)) with 

AS 1} (t) = V(M*)) 2 + (M*)) 2 + fo(*)) 2 - (41) 

The positive eigenvalue-process (41) gives the numerical realization of r(i),0 < t < 1; the 
3-dimensional Bessel bridge. Actually the sample of path of r(t) shown in Fig.l given 
in Introduction was numerically drawn by this method. In order to check the validity 
of the present numerical method to simulate the Bessel bridge, we have generated 1000 
samples of paths and studied the distribution of their maximum values = max <t<i r(t) 
numerically. The obtained result is plotted in Fig.4. There the exact curve obtained from 
Eq.(2) is also shown. The coincidence is excellent. 

B. Means and variances 

Note that the 2x2 matrix (40) can be considered as the special case of (19) with N — 1. 
If we prepare ten independent Brownian bridges, {^(t)}]^!, numerically, we can simulate 



19 



TV 


(AO 


v {N) 

N 






KIK m<J° 


KIK vW 


1 


1.251 


0.0774 


1.2533 


0.0737 


1.253314 


0.074138 


2 


1.819 


0.0732 


1.8222 


0.0746 


1.822625 


0.073194 


3 


2.262 


0.0704 


2.2677 


0.0720 






4 


2.641 


0.0664 


2.6460 


0.0692 






5 


2.979 


0.0640 


2.9805 


0.0656 






6 


3.280 


0.0624 










7 

8 


3.558 
3.817 


0.0600 
0.0556 










9 


4.057 


0.0570 










10 


4.291 


0.0543 










20 


6.146 


0.0409 










30 


7.597 


0.0396 










40 


8.790 


0.0384 










50 


9.841 


0.0364 










60 


10.806 


0.0354 










70 


11.678 


0.0310 











TABLE I: Comparison of the values of m N and v N obtained by the present numerical method 
and the exact values. The values with 'F' are read from Table 1 in [31] and those with 'KIK' are 
from Table 1 in [13]. 



the 4x4 matrix- valued Brownian bridge in 7i c (4), 

/ h(t) ^( b2 (t) + ib 3 (t)) h(t) + Mt) ±(b 7 (t) + ib 8 (t)) \ 

h(t) - Mt) ±(b 7 (t) - ib 8 (t)) -h -j= 2 (Ut) - M*)) 

\-t=(b 7 (t)-ib 8 (t)) b 9 (t) - ib 10 (t) -±(b 2 (t)+ib 3 (t)) -h(t) J 

(2) 

which is the N = 2 case of (19). By tracing the two positive eigenvalues A+ (i) = 
(\f\t),\ { 2 2 \t)),0 < \f\t) < \ { 2 2 \t),0 < t < 1, we can simulate the noncolliding paths 
of two Bessel bridges, r^ 2 )(t) = {rf\t),rf\t)) and statistical data of = max 0<t<1 rf\t) 
can be obtained. In general, we can simulate the TV-noncolliding Bessel bridges r^ N \t),0 < 
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FIG. 5: l/TV-plot of the estimated values of vq 

t < 1, by using numerical data of independently generated N(2N + 1) Brownian bridges 
{bj(t)}JSl N+1 ^ and by tracing the iV positive eigenvalues of the 2iV x 2N matrix. 
From now on, we use the notations 

»<"> = var(i/f> ) = <(<> - mf>) 2 >, 1 < * < N (42) 

for the means and variances of maximum values of paths. Table I shows the numerical 
results for the outermost paths k = N up to iV = 70, where averages have been calculated 
over 1000 samples. The present numerical results are consistent with the exact values, which 
can be read from the previous papers by Feierl [31] and by the present authors [13]. 



C. N — ► oo asymptotics of the outermost paths 

Now we study the N — > oo asymptotics of . We assume the form 

mJJ° = CoN»°, N > 1. (43) 

We have prepared a set of successive five results of the computer simulation, 
(m^\ m^^) and fitted them to the relation (43) to estimate the ex- 

ponent u and the coefficient c by the least-square fitting. We have observed that the 
estimated values of u and Co change rather systematically as the values Nj,l < j < 5, in- 
crease. Figures 5 and 6 show the dependence of the estimated values of u and c on N = N 3 . 

21 



1.45 r 



o 
O 




1.2 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 
1 //V 

FIG. 6: l/TV-plot of the estimated values of cq 




A/ 

FIG. 7: Log-log plot of vs. N for iV-noncolliding Bessel bridges. 

In these 1/iV-plots [58] of the estimated values, we made linear-fitting of the largest three 
plots as shown by the lines in the figures and obtained the values u = 0.493 and c = 1.42. 

They are consistent with Eq.(7); u = 1/2 and c = \pl = 1.414 In Fig. 6 we can see 

that the plots in the intermediate region 0.1 < 1/N < 0.2 give the values 1.25 < c < 1.3. 
They may correspond to the value y/1.67 ~ 1.29 found in the estimate (6) by Bonichon 
and Mosbah [46], which was claimed by Schehr et al. as the pre-asymptotic behavior [32]. 
Figure 7 shows the log-log plot of raffl versus N, where the curve is obtained by fitting the 
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FIG. 8: Log-log plot of vffl vs. N for iV-noncolliding Bessel bridges. 



function 



m 



(AT) 
N 



/ 2N + c 1 N~ Ul with V\ 



6 



(44) 



to the data. The fitting parameter c\ is determined as C\ ~ 0.253. 

We have also studied the asymptotics of variance v^P in N — > oo. Figure 8 shows the 
log-log plot. By the least-squares fitting, we have obtained 



~ 0.09iV- a23 . 



(45) 



The negative exponent implies i> 



(AT) 
N 



in JV -> oo. 



D. On inner paths 



Using the numerical data of paths of iV-noncolliding Bessel bridges generated by the 
method mentioned in Sec.IV-B, we can examine the statistics of the maximum values at- 
tained in the time interval (0, 1) of not only the outermost paths but also all inner paths. 
For example, the means and the variances v[ N ^ for all paths 1 < k < N are given in 
Table II for N = 10. 

Now we report our observation of the iV-dependence for the maximum values of inner 
paths. Figures 9 and 10 show the A^-dependence of the means and variances for the maximum 
values of the innermost path, h[ n \ By the log- log plots, we have obtained the following 
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k 


(10) 

mi 

K 


(10) 

vi ' 

K 


1 


0.547 


00528 


2 


0.891 


00778 


3 


1.240 


00971 

\J . \J\J U \ -L 


4 


1.598 


0.0126 


5 


1.958 


0.0147 


6 


2.337 


0.0170 


7 


2.735 


0.0209 


8 


3.168 


0.0262 


9 


3.659 


0.0320 


10 


4.291 


0.0553 



TABLE II: Numerical values of m k and v k evaluated by computer simulations 

10\ 




N 

FIG. 9: Log-log plot of m[ N ^ vs. iV for A^-noncolliding Bessel bridges. 

power-law behavior of the N- dependence. 

m[ N ^N-°- 38 , v^-N- 1 - 17 . (46) 

Figure 11 shows the dependence of on k, which is the index of paths counting from 
inner to outer. We can see that all plotted curves in Fig. 11 have a common feature as a 
function of k. This fact is clarified by Fig. 12, in which the data collapse [5] is shown by 
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FIG. 10: Log-log plot of v[ N ^ vs. N for iV-noncolliding Bessel bridges. 
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FIG. 11: Dependence of on k shown for various N. 

plotting the quantity mj^/m^ against k/N. The result implies that there is a universal 
function f(x) such that the following relation holds for sufficiently large N, 



m 



(JV) 



rn 



(JV) 
AT 



./ 



(47) 



The scaling function f(x) has two regions separated by a crossover point x*. For small 
the scaling function /(x) behaves as a linear function. On the other hand, for large 
x 3> x*, f(x) does not behave as a simple linear function. We have estimated the following 
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FIG. 12: Data collapse for the maximum values of inner particles for iV-noncolliding Bessel bridges, 
by numerical fitting. 

{x (x <ti x*), 
(48) 
a x (x ^> x*) with a ~ 5.1. 

V. CONCLUDING REMARKS 



In this paper we have reported the exact and the numerical results on the maximum-value 
distributions of paths in the A-noncolliding Bessel bridges. We have shown that the present 
maximum-value problem for a version of vicious walk model of statistical physics is related 
with the random matrix theory, the representation theory of symmetry, and the number 
theory. There are a lot of open problems. We will list up here some of them. 

In the expression (21) of the transition probability density p^ hs of the absorbing BM in 
an interval (0, h), a variable n e N = {1, 2, 3, . . . } was introduced to indicate modes in the 
Fourier expansion. When we consider the A-path systems, a set of N discrete variables n = 
(ni, n,2, . • • , %) G is introduced. Though the variables n are auxiliary, since the physical 
quantities are given by the summations over n's as shown in (21), (27), (30) and (36), we 
have seen in the derivation of Lemma 1 given in Appendix A the discrete variables n behave 
as duals of the continuous variables x = (xi,X2, ■ ■ ■ ,xn) and y = (y±, y 2 , . . . , yjv), which 
are physical variables indicating the positions of particles. The correspondence between the 
probability density of paths (14) given in the form of that of eigenvalues of random matrices 
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in class C and the distribution function of the maximum value (27) given in the form of the 
"partition function" of discrete variables implies some duality relation. The maximum- and 
minimum-value problems of watermelons without wall recently studied by Feierl [57] and by 
Schehr et al. [32] are very interesting. Systematic study will be desired to clarify the duality 
between the noncolliding path problems [18] and their extreme- value problems not only for 
bridges (i.e., excursions, watermelon configurations) [44, 45] but also for meanders (i.e. star 
configurations) [43]. 

As shown in Introduction, the distribution of the maximum value of the Bessel bridge 
= max <t<i r(t), and that of the position of Bessel bridge r(l/2) at time t = 1/2 are 
quite different from each other, because of large fluctuation of path. For example, the mean 
value of the maximum (H[^) = a/tt/2 ~ 1.2533 is much bigger than (r(l/2)) = y/2/ir ~ 
0.7979. As shown by Fig. 8, however, the variance will vanish in iV — > oo as Eq.(45). 
Then we expect 

(H$ } ) ~ <r#°(l/2)) in N - oo. 

It is known that 

(rjf \t)} ~ 2y/2Nt(l-t), < t < 1 in N —> oo 

(see, for example, Eq.(2.26) in [45], in which we should put a = and use the relation ^/x = 
lim Jv ^ 00 (r5f ) (t))/v / 2lV.) Then Schehr et al. concluded Eq. (7). By taking the proper scaling 
limit associated with iV — > oo, the fluctuation of the outermost path r^\t),0 < t < 1, 
defines the Airy process [44, 59-61], which obeys the Tracy- Widom distribution [62] at each 
time < t < 1. The present study suggests us to study the maximum value distribution of 
the Airy process. 

We have used the exact expressions (27), (30) and (36) for distributions of H$ ] given 
in Sec. Ill to check the validity of numerical calculations by our computer programs. As 
mentioned in Remark 1, the expression (27) seems to be a "partition function" of some 
discrete statistical-mechanical model with the Boltzmann weight exp(— 7r 2 |n| 2 /(2/i 2 )). The 
expression (30) reminds us the bidirectional Wronskian solutions of nonlinear equations. And 
as demonstrated below Theorem 4, the expression (36) is useful to reproduce the previous 
results reported in [12, 13, 29] for small N. Deeper understanding of these expressions in 
physics is desired. At the present stage it is not obvious how to discuss the asymptotics of 
these expressions in N — > oo. We hope that our determinantal expressions (30) and (36) will 
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be useful, since determinantal formulas play important roles in the random matrix theory 
in analyzing the large-matrix limit [65], when they are equipped with proper mathematical 
tools, e.g. orthogonal polynomials, Fredholm determinants and so on [18, 40]. 

It will be challenging open problem to analyze the maximum-value distributions for the 
inner paths in the present systems. 
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APPENDIX A: DERIVATION OF THE ASYMPTOTICS (25) AND (26) 

By the multilinearity of determinant 



{A °(i,y|x) 



det . 

l<j,k<N 



1 p -(|x| 

(2ir) N / 2 




l<j,k<N 



Here 



det [e % 



e -yj x k 



l<j,k<N l 



N 



N 





0<mi<?ri2<---<mjv j=l 




det [y 2mk ] det [x 2 

:j,k<N L J J l<j,k<N L J 



,2m fc 



28 



Now we change the variables in summation from m 3 - to fij by = rriN-j+i—N+j, 1 < j < N, 
and introduce the Schur function [48, 63, 64] 



det 

l<j,k<N 



u k +N-k 



det [xf- k ] ' 

l<j,k<N L J J 



where the denominator is the Vandermonde determinant 



(Al) 



det [xf~ k ] = TT ( Xj -x k ). 

Kj,k<N 1 3 J J.J. v J 

~ l<3,k<N 

The Schur function expansion is readily performed as (see [18, 39]) 

^)(l,y|x) = Q^-(NW)/ 2 ^ ( ^.) JJ {^-xl)(y]-yl)} 



(A2) 
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where £(fi) denote the number of parts of the partition \x (that is, the number of nonzero 
VjA < j < N). Since s„(0) = unless n = = (0, • • •, 0) e N^, and s (0) = 1, the 
asymptotics (25) is obtained. 
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where N = {1, 2, • • • } and Sn is the set of all permutations of A" items {1, 2, ... , A^}. Here 
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Therefore 
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This gives (26). 
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APPENDIX B: PROOF OF PROPOSITION 2 

From (27), using (A2) 
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where we have set k = cr(j). Since the summation XlneN^ * s taken, the above is equal to 
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where t = po a 1 and the relation sgn(cx)sgn(p) = sgn(r) was used. Then the first equality 
of (30) is proved. It is easy to confirm the second equality. 
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